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The VLT LBG Redshift Survey - IV. Gas and galaxies at 
z ^ 3 in observations and simulations 
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>-< . ABSTRACT 

CIh, We use a combination of observations and simulations to study the relationship 

JL ■ between star-forming galaxies and the intergalactic medium at z w 3. The observed 

5_( ' star-forming galaxy sample is based on spectroscopic redshift data taken from a com- 

"^ , bination of the VLT LBG Redshift Survey (VLRS) data and Keck LRIS observations 

^ ' in fields centred on bright background QSOs, whilst the simulation data is taken from 

the GIMIC simulation. In the simulations, we find that the dominant peculiar veloci- 
ties are in the form of large-scale coherent motions of gas and galaxies. Gravitational 
infall of galaxies towards one another is also seen, consistent with expectations from 
linear theory. At smaller scales, the rms peculiar velocities in the simulation overpre- 
qT ■ diet the difference between the simulated real- and z-space galaxy correlation func- 

ly^ ' tions. Peculiar velocity pairs with separations smaller than 1 h~^Mpc have a smaller 

^^ , dispersion and explain the z-space correlation function better. The Lya auto- and 

cross-correlation functions in the GIMIC simulations, appear to show infall smaller 
^U I than implied by the expected PLya ~ 1-3 (McDonald et al.). There is a possibility 

pf-j . that the reduced infall may be due to the galaxy wide outfiows implemented in the 

simulations. 

The main challenge in comparing these simulated results with the observed Keck -f 
VLRS correlation functions comes from the presence of velocity errors for the observed 
LBGs which dominate at < 1/i^^Mpc scales. When these are taken into account both 
^ , real- and z-space LBG correlation functions prefer the high amplitude of clustering 

C^ • shown by higher mass {M^, > lO^Afo) galaxies in the simulation. The simulated cross- 

correlation function shows similar neutral gas densities around galaxies as are seen in 
the observations. The simulated and observed Lya z-space autocorrelation functions 
again agree better with each other than with the /Slkq ~ 1-3 infall model. Finally, 
little direct evidence is seen in either simulations or observations for high transmission 
near galaxies due to feedback in agreement with the results of Rakic et al. However, 
these authors claim evidence for random peculiar velocities between gas and galaxies 
at small scales and strong infall at larger scales which are much less evident in our 
GIMIC simulations or indeed in the observed Keck-|-VLRS data. 

Key words: galaxies: high— redshift , intergalactic medium 



1 INTRODUCTION mation and evolution. Semi-analytical cosmological models 

require some injection of feedback to regulate the star forma- 
Tlie effect of feedback via supernovae and AGN driven winds ,. ^. .^ . i ^ i- -^ i ii ^i i ri 

tion activity m order t o limit both the number oi low-mass 

IS thought to be a key factor m the process oi galaxy tor- ,1-1 1 . JTTTTT — — 7^^ |,„„J ITTTTT — 2 — 7^ n 

*= ■' f b J and h igh-mass galaxies (jWhite fc ReeJl978l . lWhite fc FrenkI 

[l99l|). In addition, cosmologica l smoothed particle hydrody - 

namics (SPH) simulations from lSpringel fc Hernguistl (12003) 

* E-mail: pimpunyawat.tummuangpak@durham.ac.uk have shown that galactic outflow from supernova feed- 
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back is fundamental to recreating the cosmic star-formation 
history. It is also evident that simulations lacking some 
sort of feedback struggle to reproduce realistic disk galax- 
ies |Scannapieco et alj l2012l : ISchave et all I2OI0I : IWeil et akl 
Il998r i and that powerful galactic winds are required in or- 
der to produce the observed metal enrichm e nt of the IGM 
llCen fc Ostrikeij Il999l: Ixheuns et all I2OO2I : lAguirre et all 
I2OO5I : lOppenheimer fc Davdl2006l ). 

In t erms of observing the e ffects of feedback at high 
redshift, lAdelberger et alJ (|2003l . A03 hereafter) presented 
the cross-correlation between z ~ 3 galaxies and the IGM 
(as traced by QSO sightlines) and claimed an observed lack 
of absorbing gas within ~ 0.5 /i~^Mpc. They interpreted 
this as evidence of strong galactic winds removing Hi gas 
from the vicinity of these star-forming galaxies. The work 
was based on the Keck HIRES {R ~ 40, 000) spectra of 
8 background QSOs at z ~ 3 combined wit h 431 Lyman 
Break Galaxies (LBGs) from th e survev of pteidel et al.l 
(l2003h. Following the results of lAdelberger et all (|2003h . 
lAdelberger et al.l (|2005l . A05 hereafter) updated the result 
with greater numbers of galaxies, this time centred at z ~ 2. 
Based on this new sample, A05 found increase in Lya ab- 
sorption down to scales of r ~ 0.5 h~^ comoving Mpc of LBG 
positions, with no evidence for Hi gas havin g been removed 
from the vicinity of these galaxies. Indeed, ICrighton et al.l 
(|201ll ) surmised that the cross-correlation at such small 
scales would likely be affected by uncertainties in the galaxy 
redshifts in the A03 data. It is therefore still unclear to what 
extent galactic winds have an effect on the galaxy surround- 
ings. 

In addition to the above evidence for gas outflows, 
gas inflows or infall down to galaxy scales is also expected 
in simple models of galaxy formation. Gas inflow is ex- 
pected to be coherent down to the virial radius of a mas- 
sive galaxy (~ 50 kpc), below which scale the situation is 
more complicated due to shocks and the gas pressure be- 
coming more important. Gas flow infall into galaxies along 
filaments is also expected in secular models of galaxy for- 
mation where the gas accretion rate may not b e simply dic- 
tated by merging rates in a hierarchical model (JDekel et al.l 
I2OO9I ). iRakic et al.l (|2012l ) presented a study of the galaxy- 
IGM cross-correlation at ^ « 2.4 using 15 fields of the Keck 
Baryonic Survey. They saw fingers-of-god on sub- 100 kpc 
scales and evidence for infall on larger scales. 

In order to constrain simulations of galaxy formation, it 
is imperative to provide extensive observations of the IGM 
via hydrogen and metal absorption lines and thus identify 
and probe the infall and outfiow processes. As such, we are 
undertaking a large galaxy survey centred on distant bright 
QSOs in the f orm o f the VLT LBG Redshift Survey (VLRS). 
iBielbv et al.l (|201ll ) presented the first stage of the galaxy 
survey, comprising ~ 1,000 2: ~ 3 galaxi es within ~ 30' of 
2 > 3 QSOs. Using this sample, iCrighto n et al.l (|2011f ) per- 
formed a cross correlation analysis between the galaxy posi- 
tions and the Lya forest of the available QSO spectra in the 
fields, finding increased absorption within ~ 5 /i~^Mpc of 
galaxy positions. This result was consistent with the results 
of A03 and A05, but lacked the galaxy numbers to probe the 
~ 0.5 ft^^Mpc scales at which A03 claimed to see the effects 
of galaxy winds. Since then, the VLRS has been extended 
to incorporate ~ 2000 LB Gs within 9 separa te fields con- 
taining bright 2 > 3 QSOs JBielbv et al.ll2013l ). comparable 



in number t o the the only other equivalent surve ys at this 
redshift fe.g. lRakic et al.ll2012l : lRudie et al.ll2012l'). 

In this paper, we update the wor k of ICrighton et al.l 
120111 ). adding the ga l axy r edshifts of lBielbv et alj l|201ir ) 
and also ISteidel et al.l (|2003l ) in conjunction with the avail- 
able high-resolution QSO spectra in these survey fields. This 
work thus combine s the higher galaxy sampling rate of the 
ISteidel et al] (|2003l ') survey with the wide fields of the VLRS 
and provide a galaxy sample that can probe the full range 
of scales from a few hundreds of kpc to tens of Mpc. This 
large range of scales is imperative for distinguishing between 
models of gas infiow and outfiow in 2-D galaxy-gas cross- 
correlation analysis. 

Additionally, we use the Gala xies-Intergalactic Medium 
Interaction Calculation (GIMIC, ICrain et al.ll2009l ) simula- 
tion to help interpret the above correlation function results. 
The simulations are used to create synthetic Lya spectra 
and galaxies. We study both galaxy clustering and the rela- 
tionship between gas and galaxies via the (cross-) correlation 
functions ^(s) and S,{(j, tt). 

This paper is organised as follows. Observational data 
from the VLT and Keck surveys are described in section 2. 
Section 3 describes the GIMIC simulations. The simulated 
galaxy clustering results and their interpretation are shown 
in section 4, while the galaxy-IGM cross-clustering and Lya 
auto-correlation results are similarly described in section 5. 
Discussion and conclusions are in section 6 and 7 respec- 
tively. 

Throughout this work, we use a cosmology with Q.m = 
0.30, r^A = 0.70 and ho = 100 h kms"^Mpc"\ Unless stated 
otherwise all distances are comoving in units of h~^ Mpc, 
and magnitudes use the AB system. 



2 OBSERVATIONS 

The data used here are a combination of spectroscopically 
identified 2 ~ 3 star-forming galaxies and high-resolution 
spectral observations of the Lya forest of 2 > 3 QSOs. 

The gala xy data are a comb inatio n of the VLRS data 
presented by iBielbv et al] (|201ll 'l and iBielbv et all (|2013t l 
and the pub l icly a vailable Keck LBG data presented by 
ISteidel et al] l|2003h . We note at this point that the VLRS 
data offer coverage across large fields of view, effectively 
tracing the large scale structure across scales of ~ 30 — 60 
comoving Mpc tangential to the line of sight. In contrast, 
the Keck sample covers relatively small areas (~ 8 — 10 co- 
moving Mpc) , but offer higher sampling rates of the galaxy 
population. The QSO spectra with which we trace the dis- 
tribution of Hi within the fields have all been obtained from 
archival VLT UVES and Keck HIRES observations. In this 
section, we give details of all the data used and the reduction 
processes used for the QSO spectra. 



2.1 LBG Observations 

The VLRS currently provides ~ 2, 000 spectro scopic galaxy 
redsh i fts wi thin 9 fields centred on 2 > 3 QSOs (|Bielbv et al.l 
l201ll . I2OI3I '). The redshifts were obtained using the VLT 
VIMOS instrument with the LRJBlue grism giving a res- 
olution of _R ~ 180 and velocity accuracies of 5v ~ 
350 kms^^. In total the survey covers an area of ~ 2.6 
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deg^ and provides galaxy data in the foreground of the 
following 9 high redshift QSOs: Q0042-2627 (2 = 3.29), 
J0124+0044 (2 = 3.84), Q0302-0035 {z = 3.23), HE0940- 
1050 {z = 3.05), J1201+0116 {z = 3.23), PKS2126-158 
(z = 3.28), Q2231+0015 {z = 3.02), Q2348-011 (2 = 3.02) 
and Q2359+0653 (2 = 3.23). The spectroscopic galaxy sam- 
ple is predominantly limited to _R < 25 although a number 
of fainter galaxies (i? < 25.5) are present in the sample 
where slit allocation during the VIMOS observations could 
be optimised by their inclusion. 

The LBG redshifts were identified using Lya emission 
lines and interstellar medium (ISM) absorption lines where 
visible. For both the Lya and ISM features it is necessary 
to correct the ascertained redshift for intrinsic velocity ef- 
fects , due to these featu res originating in outflowing gas (e.g. 
A03. ISteidel et al.ll201ol ). As such the VLRS galaxy redshifts 
have been co r rected according to the prescription given by 
ISteidel et all (|2010h . 

The Keck survey provides a sample of ~ 940 LBGs ob- 
served using the Keck LRIS instrument. The QSOs from 
six (Q0201-fll20, Q0256-0000, Q0302-003, Q0933-f2854, 
Q2233-f 1341 and Q1422+2309) out of the 17 Keck fields 
are available to us through the public archive and taking 
only those galaxies in fields around these 6 Keck QSOs, the 
numbers of LBGs are reduced to 308. The Keck LBGs are 
limited toTZ = 25.5. 



2.2 QSO data 

We have analysed publicly available archival spectroscopy 
for 17 QSOs in the redshift range 2.5 < 2 < 4.0. The data 
are all high resolution, high signal-to-noise spectra observed 
using either the UVES instrument on the VLT or the HIRES 
instrument on the Keck telescope. Full details of the reduc- 
tion of UVES and HIRES QSO spectra for 11 of the QSOs 
used here are provided by ICrighton et al.l (|201l|) . The re- 
maining 6 spectra were all observed with the Keck HIRES in- 
strument and were reduced following an identical me t hod t o 
that used for the two Keck QSOs of ICrighton et al.l (|201ll ). 
using the makee packagqj. Briefly, this encompassed ba- 
sic flat-fielding and bias subtraction, followed by the use of 
SPIM2 to splice the echelle orders and combine individual 
observations. This involved producing template spectra con- 
structed by combining the individual observations, masking 
bad regions of the CCDs and rescaling. A template was ap- 
plied to rescale the original observations. We divide out the 
continuum for each individual observation, then multiply 
this normalised flux by a continuum fit to the template. Af- 
ter scaling each order of each observation individually, we 
combined them to get the final spectrum. 

In addition to the publicly available QSO spectra, we 
also include a spectrum from our own observations using 
the X-Shooter instrument on the VLT (Q2359-I-0653). This 
data was reduced using the X-Shooter pipeline package, see 
Bielby et al. (in prep) for details. The full list of QSOs used 
in this study is provided in Table [l] 



http : //www2 . keck, hawaii .edu/inst/hires/makee www/ 



3 GIMIC SIMULATIONS 



3.1 Overview 



We use gas dynamical simulations to create both Lya spec- 
tra and galaxies to compare with the observational data. Our 
main aims are to study the real and redshift-space correla- 
tion functions to see if we can detect the effects of peculiar 
velocities to understand more about gas outfiow and infall 
into galaxies, for (a) LBG-LBG pairs (b) the LBG-Lya for- 
est and (c) Lya-Lya pairs. The results will then be used to 
interpret the observable 1-D and 2-D correlation functions 
^(s) and ^{o, tt) in terms of both simulation and observa- 
tional results. 

For this purpose we use the GIMIC simulation, which is 
a cosmological hydrodynamical re-si mulation of sele cted vol- 
umes of the Millennium simulation (|Springelll2005l '). GIMIC 
is designed to overcome the issues in simulating large cos- 
mological volumes (L > IQQh^^ Mpc) at high resolution 
(rrigas ^ 10' h~^ Mq) to 2 = by taking a number of smaller 
regions with 'zoomed' initial conditio ns (jFrenk et al.lll996l : 
iPower et al]|2003l : iNavarro et al.ll2004l ') . These individual re- 
gions each have approximate radii of 18/1"^^ Mpc outside of 
which the remainder of the Millennium simulation volume 
is modelled with coUisionless particles at much lower reso- 
lution. 

In the work presented here we use only the 'Oa' GIMIC 
region, which is identified as having a mean density at 
2 = 1.5 equal to the mean density of the simulation as a 
whole at that epoch. In addition, we use only one snapshot 
of this region, chosen to be at a redshift of 2 = 3.1 in order 
to provide a suitable comparison to our 2 ~ 3 observed pop- 
ulation of star-forming galaxies. GIMIC was run using the 
TreePM SPH co de GADGET3 , which is an update of the 
GADGET2 code (|Springelll2005h . The cosmological parame- 
ters adopted were: i},n = 0.25, Q\ — 0.75, fit, = 0.045, ho = 
100 h kms"^Mpc"\ h = 0.73, ag = 0.9 and n^ = 1 (where 
Us is the spectral index of the primordial power spectrum) . 
We note also that the GIMIC simulation does not incor- 
porate the effects of AGN feedback but, importantly, does 
implement a prescription for supernovae driven feedback. 

Although the GIMIC volume is small, it still can be use- 
ful for studying 2 — space distortions. This is because in the 
simulations we have knowledge about the real-space corre- 
lation functions and so we do not need an ensemble average 
over a large volume to guarantee spherical symmetry. 



3.2 Simulated galaxy population 

The galaxy population is identified in the simulation based 
on first identify ing the dark mat ter halos using a Friends 
of Friends (FoF. IPavis et al.lll985l ') algorithm. A group find- 
ing algorithm then locates the nearest dark matter halo for 
each baryonic (gas or star) particle and ide ntifies the par- 
ticle with this halo. Th e SUBFIND algorithm (jSpringel et al.l 
I2OOII : iDolag et aLll2009l ') is then used to identify self-bound 
sub-structures within the halos, to which star particles are 
associated and defined as galaxies. 

Given a limiting radius for our galaxy sample within 
the GIMIC volume of 16 h~^ Mpc radius simulation vol- 
ume and a stellar mass limit of M* ^ 10* M0 we obtain 
a total of 4,070 galaxies from the snapshot at 2 = 3.06 in 
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QSO 


R.A. 


Dec. 


z 


Mag 


Instrument 




J2000 








Q2359+0653 


00:01:40.6 


+07:09:54 


3.23 


V = 18.5 


X-Shooter 


Q0042-267 


00:44:33.9 


-26:11:19 


3.289 


B = 18.5 


HIRES 


J0124+0044 


01:24:03.8 


-1-00:44:32 


3.83 


9 = 19-2 


UVES 


Q212904.90-160249.0 


21:29:04.90 


-16:02:49.0 


2.90 


19.2 


HIRES 


[iy//O91]0043-265 


00:45:30.470 


-26:17:09.20 


3.44 


18.34 


HIRES 


Q0201+1120 


02:03:46.7 


-1-11:34:45 


3.610 


20.1 


HIRES 


Q0256-0000 


02:59:05.6 


-1-00:11:22 


3.364 


18.2 


HIRES 


LBQS 0301-0035 


03:03:41.0 


-00:23:22 


3.230 


17.6 


HIRES 


Q0302-0019 


03:04:49.9 


-00:08:13 


3.281 


17.8 


HIRES 


QSO B0933+289 


09:33:37.2 


-1-28:45:32 


3.428 


17.8 


HIRES 


HE0940-1050 


09:42:53.5 


-11:04:25 


3.06 


B = 17.2 


UVES 


J1201+0116 


12:01:44.4 


-1-01:16:11 


3.233 


9 = 17.7 


HIRES 


PKS2126-158 


21:29:12.2 


-15:38:40 


3.268 


V = 17.3 


UVES 


Q2231-0015 


22:34:08.9 


-1-00:00:01 


3.02 


r = 17.3 


UVES 


Q2233+136 


22:36:27.2 


-1-13:57:13 


3.209 


18.7 


HIRES 


Q2348-011 


23:50:57.9 


-00:52:10 


3.023 


r = 18.7 


UVES 


Q1422-I-2359 


14:24:38.1 


-1-22:56:01 


3.620 


16.5 


HIRES 



Table 1. List of QSOs used in this study. 
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Figure 1. The distribution of A/* > 10^ M© simulated galaxies 
in real-space (pink asterisks) and redshift-spacc (diamonds). 
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Figure 2. The histograms of pairwrise velocities (Aiiz) (pink his- 
tograms) from the full sample and from the s < 1 h~^ Mpc 
sample of galaxy pairs (blue dashed histograms) for the high and 
low mass galaxy samples. 



the Oct density region. A number of studies have investi- 
gated the mass distribution of LBGs at z « 2 — 3, suggest- 
ing median stellar masses for galaxies selected in this way of 
M^ ~ 1.4x10^° Mq with a scatter of ~ logio(M/M0) ~ 0.5 
(jShaplev et al.l 12009 ). However, at such masses, the galaxy 
population in GIMIC is too small for the analysis we wish to 
perform and therefore a compromise is required. We there- 
fore choose a mass cut of M* ^ 10^ M0 to provide our 
closest feasible estimate of the LBG population. The effect 
of choosing a relatively low mass cut is somewhat mitigated 
by the relatively high value of cts used in the GIMIC sim- 
ulations, which acts to increase the simulated clustering. 



Thus, although the mass cut is a little low, the clustering of 
the simulated sample should be comparable to the observed 
LBG population, given the large erg in GIMIC. Doing so 
gives a sample of M* > 10^ M© galaxies with r ^ 16/i^^ 
Mpc of the simulation volume, giving a space density of 
Pg ^ 5 X 10~^ cMpc~^ and 287 galaxies. This is marginally 
higher than, but still com parable to, the observe d number 
densities of LBG galaxies (jAdelberger et al.ll2005l ). 

The distribution of M* ^ 10^ M© galaxies in real- (di- 
amonds) and redshift-space (asterisks) is shown in Fig. [l] 
Throughout this paper, we will assume that the x and y co- 
ordinates give positions tangential to the line of sight, whilst 
z coordinates give the line of sight positions, either in real 
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Figure 3. The projected position of A/* ^ 10 Mq galaxies (di- 
amonds) and 200 Lyo sightlines (circles) on the XY plane. 



or redshift-space. Fig. [T] illustrates the measured positional 
shifts in the z-axis given by the peculiar velocities of the 
galaxies within the simulation. It is evident from this plot 
that there is a consistent 'bulk' motion directed in the pos- 
itive z-direction. Measuring the distribution of the galaxy 
velocities, we find an average velocity {vz) = 93 kms~^ with 
a standard deviation of 128 kms~^ for the M* ^ 10^ M© 
galaxy sample and {v,) — 4-94 kms~^ with a standard de- 
viation of 125 kms"^ for the A/* ^ 10* M© galaxy sample. 
We show the pairwise velocity (Aii^) distributions (pink 
histograms) of galaxies in Fig. [2] For the IS/'h ^ lO'' M© and 
A/* ^ 10 M0 galaxy samples we find a standard devia- 
tion of (Jv — 177 kms~^ and a„ = 173 kms~^ respectively. 
The blue dashed histograms show the distribution for only 
those pairs within r = 1 h~ Mpc of each other, thus iso- 
lating the intra-halo velocity dispersion and excluding the 
effect of the halo-halo velocity dispersion. This is impor- 
tant when considering the effect of the velocity dispersion 
on the galaxy-galaxy clustering measurement. The standard 
deviations of the pairwise velocities for pairs within s ^ 1 
/i-^Mpc are 134 kms^^ and 100 kms"^ for Af* ^ 10** Mq 
and Mi, J5 10* Mq galaxies respectively. None of these stan- 
dard deviations include redshift uncertainties due to mea- 
surement errors that affect the observed galaxy redshifts. 

3.3 Simulating Lya forest spectra 

We next generated spectra along the 2-direction through the 
GIMIC volume to compare with the observations. The sight- 
lines were extracted using SPECWIZARe|3. The transmission 
is given by, T = e^^, where r is the optical depth along 

^ Developed by J. Schaye, C. Booth and T. Theuns, refer to 
iTheuns et al] lll998l) for details 



Figure 5. The distribution of LOS mass weighted peculiar veloc- 
ities measured within each pixel in each of the GIMIC simulated 
spectra. This illustrates the underlying dynamics present in the 
spectra. The LOS peculiar velocity distribution shows a mean pe- 
culiar velocity of [v^.) = 110 kms~^ with a standard deviation of 
120 kms-i. 



the line-of-sight. We use a spectral resolution FWHM of 
7.0 kms~^ to convolve each spectrum, a signal-to-noise of 50 
per pixel, and pixels of width 2.8 kms~^ which are typical 
values of our UVES and HIRES QSO spectra. The sightlines 
were generated with random X, Y positions, then parallel to 
the Z— axis. We constructed 200 sightlines, with this number 
providing an approximate optimum between maximising the 
statistics available and minimising the over-sampling of any 
region of the volume. Each sightline was constrained to not 
extend beyond 16 h~^ Mpc from the centre of the GIMIC 
volume in order to avoid any edge effects in terms of the 
gas extent. The average transmission, Tr for real space is 
0.69 while the Tz for redshift-space is 0.72. An explanation 
for this difference is that infall of saturated absorption lines 
towards each other in redshift-space may result in an overall 
increase in the measured flux. This will cause the average 
transmissivity over the full spectrum to increase in redshift- 
space as seen. Some hint of this effect can be seen in Fig. [3] 
in which we show a number of examples of the flux from 
each sightline compared in real and redshift-space. 

Using SPECWIZARD, we calculate the mass weighted 
line-of-sight (LOS) peculiar velocities for each pixel in our 
200 spectra. The distribution of the peculiar velocities is 
given in Fig. [5] As with the galaxy population, the gas 
traced by the simulated spectra shows the bulk motion in 
the positive z-direction, with a mean peculiar velocity of 
{vz) = 110 kms"'^ and a standard deviation of 120 kms~^. 
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Figure 4. Examples of absorption spectra from simulated QSO sightlines. The blue solid line shows 
line shows e~^ in redshift-space.The scale is measured in h~^ comoving Mpc. 
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4 GALAXY CLUSTERING 



4.1 1-D correlation function 



4. 1.1 Estimator 

iBielbv et alj (|2013l ') presented a clustering analysis of the 
LBG data used in this study (combining the VLRS and Keck 
data). In this sectio n, we compar e the o bserved galaxy clus- 
tering presented bv lBielbv et al.l (|2013l l to results obtained 
using the galaxy population within the GIMIC simulation. 
In so doing, we may validate how representative the GIMIC 
galaxy population is of the observed 2 « 3 LBG population 
in terms of intrinsic clustering properties and the effects of 
the galaxy velocity field on the galaxy clustering. 

We calculate the real and redshift-space functions, ^(r) 
and €,{s), of the GIMIC z — 3.06 galaxy samples using the 
estimator: 



ar) 



Nr {DD{r)) 
Ng {DR(r)) 



(1) 



where {DD{s)) is the average number of galaxy-galaxy pairs, 
{DR{s)) is the number of pairs of galaxy-randoms at the 
separation, r, where we use r to represent the separation in 
real space and s the separation in redshift-space. The factor 



Ng 



is the ratio of the number of random to data points. 



We estimate errors on the auto-correlation results using 
jack-knife estimates based on splitting the simulation into 



equal volume octants and excluding each octant in turn to 
create 8 jack-knife realisations of the data. The correlation 
functions are then fit using a power-law of the form of: 



ar) 



(2) 



where 7 is the slope of real-space clustering, (,{r), and ro is 
the real-space clustering length. 



4-1.2 Simulated real-space galaxy correlations 

Fig.[6]shows the results for the simulated galaxy-galaxy cor- 
relation function with (a) M* ^ 10^ M0 and (b) Af* > 
10* Mq simulated galaxies. The blue diamonds show re- 
sults from galaxies in redshift-space while the pink aster- 
isks show results from galaxies in real space. The integral 
constraint, X, is included in the data in order to compen- 
sate for the effect of the limited field sizes (as described in 
iBielbv et al.ll2013h . The estimated integral constraints are 
1 = 0.21 and 1 = 0.11 for M* > 10^ Mq galaxies and the 
full sample respectively. The pink lines represent power-law 
fits to the real-space correlation function based on Eqn. [2] 
with 7 = 1.56 ± 0.26 and ro = 4.16 ± 1.16 h'^Mpc in 
the M* ^ 10^ Mq galaxies' plot while 7 = 1.52 ± 0.10, 
ro = 2.41 ± 0.24 h'^Mpc in M* ^ 10* Mq galaxies' plot. 
These power-law results give good fits to the real-space clus- 
tering results and we note that there is little sign of a double 
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Figure 6. (a) Galaxy correlation functions for 287 simulated M^. ^ 10® Mq galaxies compared to the observed LBGs. Real space - 
pink asterisks, redshift-spacc - blue diamonds, LBG data - brown triangles. Also shown are fits to simulated real space (pink line) with 
7 = 1.56 and ro = 4.16 fc~^Mpc. The real-space fit to the Keck+VLRS data is shown by the grey dot-dashed line. The pink dashed 
line is the predicted ^(s) assuming the real-space fit with the full pairwise velocity dispersion of (lu^) ~ ^''"^ kms~^. The blue solid 
line is the same except with the s < 1 h~^ Mpc pairwise dispersion of (tc^) = 134 kms~^. The brown line is the same model now 
with (ui^) = 420 kms~^, with velocity errors added to allow comparison with the Keck-|-VLRS LBG £,{s). (b) The same for 4,070 
simulated Af* ^ 10* TsAq galaxies with real space fit 7 = 1.52, ro = 2.41 /i~^Mpc. The ^(s) predictions now assume the appropriate 



pairwise velocity dispersion of (tu^) 



1/2 



173 kms ^ (pink dashed line) and (to^) = 100 kms ^ (blue solid line). Bottom panels 



present £,{s)/S,{r) with jack-knife error bars. The dotted line represents the predicted Kaiser boost (Kaiser et al. 1987) with (a) /3 = 0.35 
giving C(s)/5(r) = 1.26 for M* ^ 10^ M© galaxies and the Kaiser boost with (b) /3 = 0.53 giving 5(s)/C(r) = 1.41 for M« ^ 10* Mq 
galaxies. 



power-law or two-halo break in the clustering for either of 
the samples. The contribution of the single halo term in 
galaxy clustering is generally seen to be more significant 
with higher galaxy stellar mass however and so the lack of 
any break may be due to the mass constraints used here. 
This result should clearly be checked in bigger simulations. 



4-1-3 Simulated ^{s)/^{r) and mfall 

In the lower panels we show the ratio between the real 
and redshift-space clustering results in order to highlight 
the signatures of redshift-space distortions in the redshift- 
space correlation function. Here the errors are again con- 
structed from the jack-knife realisations. At scales above 
s ~ 1.5 — 2 h~^Mpc, we see the effects of dynamical infall, 
which acts to boost the clustering sig nal in the red s hift-space 
meas urement. From linear theory (IKaiseiJ 119871 : iHamiltonl 
1 1992 ) we expect to see the 'Kaiser boost', which is given 
given by: 



C(s) = (l + |/3 + i/3^)C(r) 



(3) 



where /3 is the dynamical infall parameter. For galaxies, /3 = 
OP'^ /b where b — \/^gg/S,mm is the linear bias parameter. 
At z « 3, we proceed via the volume averaged clustering 



amplitude, ^(8), (see eqns 17, 18 of lSielbv et"aLll2013l ). For 
the Mi, ^ 10' M0 simulated galaxy case, the above power- 
law parameters fitted to ^(r) give ^g(8) = 0.75 ± 0.05 which 
with Cba/(8) = 0.088 gives bias 6 = 2.80 ± 0.18. Taking 
n™(2 = 3.0) = 0.98 gives P = 0.35 ± 0.02. So the expected 
value of the Kaiser boost is 1.26 ± 0.12 compared to the 
observed amplitude ratio at 1 < s < 8 /i^^Mpc of 1.17±0.06. 
Assuming the power-law fitted to ^(r) for the set of 
M* > 10* Mq galaxies, we find ^9(8) = 0.33 ± 0.02, giving 
6 = 1.85 ± 0.12 and (3 = 0.53 ± 0.03. The predicted Kaiser 
boost is therefore ^is)/^{r) = 1.41 ± 0.14. At separations 
of 1 < s < 8 /i~^Mpc, we find a mean amplitude ratio of 
1.26±0.03, consistent at the ~ la level with a Kaiser boost, 
although marginally lower than predicted by linear theory. 
We conclude that for the simulation, the ratio, (.{s)/£,{r), is 
in reasonable agreement with linear theory infall predictions 
for both high and low mass galaxies. 



4- 1-4 Simulated galaxy correlations and velocity dispersion 

At smaller separations (r < 1 h^^ Mpc) for both high- and 
low- mass simulated galaxies, the galaxy-galaxy ^(s) in red- 
shift space has lower clustering than £,{r). This turn-over of 
the real space correlation function is the result of z-space 
smoothing due to the peculiar velocity dispersion, (w^) 
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We model the effects of the velocity dispersion on the clus- 
tering results using a gaussian distributio n to the veloc- 



ity d i spersion, following prev ious work (e.g. iHawkins et al.l 
I2OO3I : Ida Angela et al.l |2005| ). Using the pairwise velocity 
dispersions derived from Fig. [2] for the two samples (i.e. 

{wiy^ = 177 kms"' and (w?)'''^ = 173 kms"' for the 
M* ^ 10'-* Mq and M* ^ 10* Mq samples respectively 
- pink dashed lines in both panels), we find that the re- 
duction of the real-space clustering at small scales is over- 
predicted compared to the measurements of ^(s). As illus- 
trated in Fig. [2] however, we note that the measured pair- 
wise velocity dispersion is separation dependent. The ob- 
served discrepancy is therefore likely the result of the effect 
of small scale peculiar motions on the clustering function 
being dominated by galaxies within ~ 1 h~^ Mpc of each 
other, whereas the initial pairwise velocity histogram pre- 
sented in Fig. [2] includes pairwise velocities between galax- 
ies across all separation scales within the simulation. If 
we thus limit the histogram of pairwise velocities to only 
those pairs within 1 h~^ Mpc of each other (dashed his- 
tograms in Fig. [5]), we retrieve pairwise velocity dispersions 
of {wlf^ = 134 kms"^ and {wlf^ = 100 kms"^ for the 
Mi, ^ 10^ M0 and M* > 10* M© samples respectively. Us- 
ing these values in the redshift-space distortion model, we 
find improved agreement between the model (solid blue line 
in Fig.[6ll and the galaxy auto-correlation function measured 
from the GIMIC simulations. Ultimately, the appropriate ve- 
locity dispersion for modelling the redshift-space distortion 
effects on the galaxy clustering, is the velocity dispersion 
present within groups, whilst the peculiar velocity measured 
from the simple histogram case included the imprint of the 
velocity dispersion of galaxy groups as well as the dispersion 
within groups. Taking the histogram of only pairs of galax- 
ies within ~ 1 h^^ Mpc of each other effectively measures 
the intra-group peculiar velocities. We conclude that ^(s) 
is better explained on sub-Mpc scales with the intra-group 
velocity dispersion appropriate for these scales. 



4.2 Simulated and observed correlation functions 
compared 

iBielbv et al.l (|2013l ) report the best fit scale-length and slope 
for the observed Keck -I- VLRS LBG-LBG semi-projected 
Wp{(j) for the data is ro = 3.83 ± 0.24 h^^ Mpc with a slope 
of 7 = 1.60±0.09. Within the reported errors, the clustering 
of our Mi, ^ 10^ Mq sample reproduces the observed survey 
clustering very well in terms of both clustering length and 
slope. As would be expected, the Mi, ^ 10* M© sample gives 
a somewhat lower clustering length than the observational 
data, but does at least have a consistent slope within the 
quoted errors. 

We now apply the measur ed (wz) 
servations of iBielbv et al.l (l2013h to our correl ation func- 
tions measure from GIMIC. IBielbv et al.l (|2013h measured 

{w1) = 420 kms~^, which includes both the intrinsic 
velocity dispersion and the velocity errors o n measuring 
the ga laxy redsliifts. The measured ^(s) from IBielbv et al.l 
l|2013h is shown in Fig. [6] (brown triangles) and a model 
based on the GIMIC ^(r) combined with the observational 
(w1) = 420 kms~^ is given by the brown solid line. By 
introducing the observationally measured pairwise velocity 



(w?) from the ob- 



errors to the GIMIC M* > 10^ M© result, we find that the 
GIMIC clustering measurement reproduces the measured 
LBG clustering well. 



4.3 2-D correlation function 

We now turn to the 2-D galaxy auto-correlation functions in 
order to further investigate the impact of galaxy velocities 
on clustering measurements within the simulation. In the 2- 
D correlation function, ^(ct, tt), we parameterise the line of 
sight separation between two galaxies by tt and the trans- 
verse separation by a. We calculate ^(cr, tt) using the same 
methods as used for the 1-D correlation functions and with 
the same samples. 

Fig. [7] show the 2-D galaxy-galaxy correlation function, 
i{o,-n), for Mi, > 10^ Mq and M* > 10* Mq simulated 
galaxies and the results of fitting each of these for (to^) 
and /3, using the 1-D measured clustering lengths and slopes. 
Panels (a) and (b) show the measured correlation function 
for the Mi, ^ 10^ Mq galaxy sample in real- and redshift- 
space respectively while panels (c) and (d) show the correla- 
tion function for the M* ^ 10* Mq galaxy sample. Note that 
the previously determined integral constraints are applied in 
these plots. 

Panels (e) and (f) show chi-square contours based on fit- 
ting the Mi, ^ 10® Mq galaxy ^(o", tt) for velocity dispersion 
and P, based on the real-space ^(r) (i.e. ro — 4.16 /i~^Mpc 
and 7 = 1.56). The best fit for real space was j3 = 0.02 and 
velocity dispersion — 200 kms^^ with reduced x^ = 0.7. As 
expected, there is little significant effect of peculiar velocities 
seen in real space. In redshift-space, we found /3 = 0.02 and 
the velocity dispersion = 100 kms~^ with reduced y^ = 0.7. 
Here we are expecting /3 ~ 0.35 given the bias, 6 = 2.80, and 
this is significantly different from what we find in the simu- 
lation. But the value of /3 = 0.02 from fitting ^(ct, tt) is also 
lower than the /? = 0.24±0.1 implied by the ratio of ^(s)/^(r) 
above which is in much better agreement with the expected 
13 = 0.35. For M* ^ 10* M© galaxies, we fit the model 
with 7 = 1.52 and ro = 2.41 /i"^Mpc. The best fit for real 
space gives /3 = 0.15, velocity dispersion — 80 kms^^, and 
reduced x^ = 0.4, with /? consistent with zero as expected 
and /3 = 0.45, velocity dispersion = 120 kms~^ with reduced 
X^ = 1.1 for redshift-space as shown in panels (g) and (h). In 
this case, we are expecting /3 « 0.53 given the bias, h = 1.85. 
Thus here the fit to ^(cr, tt) is in much better agreement with 
linear infall theory. Indeed, a model that assumes the above 
^(r) parameters with j3 — 0.53, (wl) = 100 kms~^ is only 
rejected by the data at the la level. We conclude that apart 
from the ^[a, tv) result for the M* ^ 10® MQgalaxies, the 1-D 
and 2-D correlation results for both high and low mass sim- 
ulated galaxies are generally in reasonably good agreement 
wit h linear infall theory . 

IBielbv et al.l (|2013 ) presented the results for ^{a, tt) for 
the VLRS-|-Keck LBGs. Again the velocity errors domi- 
nated the fits which gave j3{z = 3) = 0.38 ± 0.19 and 



w^ = 4201^60 kms~^ when ro = 3.83 ± 0.24 h'' Mpc and 
7 = 1.60 ± 0.09 was assumed for ^(r). The value of /? is in a 
range consistent with that measured above for the simulated 
galaxies. 



VLT LBG Redshift Survey IV 9 







500 






e 




500 






f 




400 










400 










A 300 

V 

200 










^A 300 

V 

200 


m', 


) 






100 










100 


f 










0.2 


0.4 0.6 


0.8 1 





0.2 


0.4 0.6 


0.8 1 









/^ga, 










Pm 







500 

400 

A 300 

s" 

V 

200 

100 





500 




400 


"a 


300 


V 


200 




100 



0.2 0.4 0.6 0.8 1.0 



0.2 0.4 0.6 0.8 1.0 

^gal 



Figure 7. The auto-correlation function ^(o-, tt) results from simulated galaxies. (a)-{b) are the results from M* ^ 10® Mq galaxies 
in real- and redshift-space respectively. (c)-{d) are the results from M* ^ 10* Mq galaxies, {e)-(f) are the pairwise velocity dispersion 
{{wl) )-infall parameter (/3) Ax^ contours for g((T, tt) with M* ^ 10® Mq galaxies, fitting to ^(ct, tt) writh s < 12 h~^'hA\>c. The best 
fit for (e) real-real gives /3 = 0.02, («)?) = 200 kms~^, assuming rp = 4.16 /i~^Mpc and 7 = 1.56. The best fit for (f ) z — 2 gives /3 = 
0.02, {wV^ = 100 kms^^. (g)-(h) show similar contour plots with (e)-(f) but with Af* ^ 10* Mq galaxies. The best fit for (g) real-real 
gives /3 = 0.15, ('U'^) =80 kms"'-, assuming ro = 2.41 h~^M\>c and 7 = 1.52 while the best fit for (h) z — 2 gives /3 = 0.45, (i«^) 
= 120 kms-i. 



5 GALAXIES AND THE IGM 

As discussed earlier, the relationship between the galaxy 
population and the IGM is key to understanding galaxy 
growth and evolution. Galaxies require large halos of gas in 
order to grow to the large masses we observe at the present 
day, whilst the supply and regulation of the flow of gas into 
galaxies dictates the distribution of galaxy masses we ob- 
serve. 

From observations of galaxy winds with speeds of 
> 300 kms^^ for the LBG population (e.g. via the off- 
sets nebulae and inter-stellar medium spectral features), 
it is evident that outflo wing materia l exist s in these star- 
form i ng galaxies (e.g. IPettini et al.l I2OOII : IShaplev et al.l 
I2OO3I : iBielbv et"aLll201ll ). A number of authors have thus 
attempted to distinguish the effects of such outflows on the 
distribution of gas around the z ~ 2 — 3 star-forming galaxy 
population via the Lya forest obser ved in the spectra of 
background sight l ines (e.g. A03, A0 5, ICrighton et al.ll201ll : 
iRudie et al.ll2012l : iRakic et al.ll2012l ). 

In this section, we therefore perform an analysis of the 
cross correlation between galaxies and the Lya forest within 
the GIMIC simulation. We apply the same dynamical mod- 
els as in the previous sections to the cross-correlation analy- 
sis. However, in the case of the galaxy-gas cross-correlation 
the relation between redshift and real-space correlations will 
become i.(s)/^{r) = (1 + \( B„ai + PLyJ + ^PgaiPLyJ (see 



This is because of the non-linear relation F = exp(— r) be- 
tween Lya transmission and optical depth, r, where most 
of the physics in the Lya forest is contained in r. Ac- 
cording to McDonald et al (2000, 2003) the infall param- 
eter Phya = ^'rn ^ bri /bs and b„ and bs have to be de- 
fined from simulations. iMcDonaldl (120031 1 found results for 
/Slj/q ~ 1 — 1.6 depending on the resolution of the simula- 
tions. We therefore took /3Ly„ = 1.3 as our estima te of the 
gas dynamical infall parameter. iMcDonaldl (|2003r ) did not 
use the z-space distortion techniques used here so this and 
the fact that we are using higher resolution SPH simulations 
makes it interesting to check whether linear theory wit h their 
I^Lya fits our simulated data. IMcDonald et al.l (|2000l ) argue 
that the form of the fiux correlation function is proportional 
to the mass corre l ation function in the linear regime. Fol- 
lowing [McDonald] (|2003l ), on this basis we shall assume that 
we can take account of 'finger-of-God' velocity dispersions 
in the usual way by convolving the transmission correlation 
function with a Gaussian of the appropriate dispersion. 



iMountrichas fc Shanksll2007l ) . The linear bias of the gas ob- 
tained from b^ = S,Lya-Lya/^m,m \s b ^ 0.3 but this is not 
the bias required to assess the effect of gas infall via /3Ly„ • 



5.1 Keck+VLRS LBG-Lya cross-correlation 

5.1.1 ID cross correlation, {T{s)) 

To measure the velocity dispersion of Hi gas that caused an 
increase in Lya absorption at the redshifts, we calculated 
the Lya fiux along QSO sightlines. We first calculate the 
transmissivity of the Lya forest by 
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T- 



f 

J con 



(4) 



where / is the observed flux and fcon is the continuum in 
flux level from Lya forest. 

Following A03, for the observational data we normalised 
the QSO transmissivities for QSOs at 2: ~ 3 by dividing the 
measured flux transmissivities, T, by 



T{z) = 0.676 - 0.220(z - 3) 



(5) 



where z is the redshift of a given pixel. This normalisation is 
used to correct for the evolution of flux transmissivity with 
redshift since the transmissivity is higher at low redshift. To 
avoid the effect of contamination of Lya absorption lines, 
we cut out the spectrum below the Ly/3 emission. Thus the 
spectrum between the Ly/? and Lya is only used in this 
calculation. We also excluded the wavelength within 20 A of 
the intrinsic Lya emission to prevent the proximity effects 
from the QSOs. 

We then use the transmissivity of the Lya forest as 
calculated above to perform the Lya-LBG cross correlation 
function. The Lya-LBG cross correlation function is calcu- 
lated from 



{T{s)) = 



{DTjs)} 
Nis) 



(6) 



where {DT{s)) is the number of galaxy-Lya pairs weighted 
by the normalised transmissivity, T/T{z), for each sepa- 
ration. N{s) is the number of LBGs that contribute to 
the cross-correlation function at each separation. This is a 
slightly different approach from A03 who use a correlation 
function based on the Landy-Szalay estimator. 

In Fig. [S] we present the LBG-Lya transmissivity 
correlation function from Keck, VLRS, and combined 
VLRS-I-A05 sample compared to the results from A03 and 
A05. The results from Keck, VLRS and combined samples 
are shown by diamonds, asterisks and filled circles respec- 
tively. The triangles and squares show the results of A03 and 
A05, respectively. At distances greater than 3 /i^^Mpc, we 
found good agreement with both A03 and A05. We note in 
passing that our own reductions of the Keck sample HIRES 
data gives a consistent result with the A03 Lya-LBG sam- 
ple. In the case of Fig. [S] we continue to use the published 
A03/A05 results with their larger number of QSOs. The 
measured {T{s)) increases with separation and reaches the 
mean value of T = 0.77. At separations below 3 h~^ Mpc, 
the combined sample seems to have the same trend as A05, 
with no evidence for a turn-up at s < 1 h~^Mpc, a fea- 
ture that was thought to be evidence for feedback. With the 
larger sample o f LBG s close to QSO sightlines compared to 
ICrighton et al.l (|201ll ). we have now strengthened the ev- 
idence against feedback strongly affecting Lya absorption 
on Mpc scales around galaxies. 



5.1.2 2-D cross correlation, ^(cr, tt) 

We now use the full VLRS sample and Keck LBG-Lya 
datasets to measure the 2-D LBG-Lya cross-correlation, 
^(ct, tt). By combining these two surveys, we can compare 
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Figure 8. The mean Lya-Galaxy transmissivity as a function of 
distance from an average z ^ 3 galaxy. Blue filled circles show 
VLRS-I-A05 result. Brown asterisks and pink diamonds show re- 
sults from VLRS and Keck respectively. A03 and A05 are also 
represented by squares and triangles. The pink dot-dashed line 
shows the GIMIC z — z result. Also shown is the fit to the 
M^, > 10* Mq GIMIC Lyo-Galaxy (T{s)) from Crighton et al. 
(2011), with so = 0.3 h-^Mpc and 7 = 1.0 (dot-dashed red line). 
The solid blue line shows this same model convolved with veloc- 
ity errors of {'w'^} = 280 kms"-"^. The bottom panel shows the 
number of LBGs contributing to each bin in each sample. 



the correlation functions in a wider range of separations than 
would otherwise be possible (the VLRS giving 2 — 3x the 
coverage in the a scale than the Keck data). The LBG-Lya 
(,{a, n) from Keck-f VLRS sample is presented in Fig. |9] 

There is little sign of any 'flnger-of-god' effect at small 
a scales (cr < 1 h~^Mpc), i.e. in the form of the clustering 
power being extended in the tt direction. We also present 
Ax^ contours from fltting LBG-Lya £,{a, tv) for velocity dis- 
persion and Pi^ya assuming 7 = 1, sq = 0.3 /i~^Mpc, and 
Pgai ~ 0.38. We assumed 7 and so from fitting LBG-Lya ^(s) 
using ^(s) = (s/so)~^. As shown in Fig. [S] this model re- 
produces the observed LBG-Lya ^(s) results well. We shall 
therefore take our starting model for the cross-correlation 
(,{r) by scaling the above ^(s) downwards by the expected 
Kaiser boost of 1.5. We checked that this factor gave the 
best reduced x^ since we are aiming to fit the shape rather 
than the amplitude of ^{a, n) here. The best fit to ^(cr, tt) was 
I^Lya = 0.02, velocity dispersion = 20 kms~^, and reduced 
X^ = 0.65. 

We shall return to these Keck-|-VLRS results to com- 
pare with the results from the GIMIC simulations described 
below. For the moment, we note that the theoretically ex- 
pected value of the infall parameter is ^Lya = 1-3 and the 
expected velocity dispersion, including the effects of LBG 
velocity error and intrinsic dispersions for gas and galax- 
ies, is (w'iy^^ = V2972 + 1202 = 320 kms"\ The observed 
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Figure 9. (a) Lya-LBG ^{cr, n) from Keck+VLRS sample, (b) The pairwise velocity dispersion - infall parameter Ax^ fitting contours 
for the Keck+VLRS data, fitting to ^(cr, tt) with s < 8 h~-'^Mpc. The best fit gives iBj^ycc = 0.02, (ui^) =20 kms"-"^, assuming sq = 
0.3 h^^Mpc, 7 = 1 (before correction for Kaiser boost), and Pgal = 0.38. 



value of Phya = 0.02 is about 3cr lower than the expected 
value (see Fig. |9] panel b). 



5.2 Lya-Galaxy ^(s) from simulations 

As with the data, we compute the galaxy-Lya cross corre- 
lation using the methods described above. We note however 
that the renormalisation to 2 = 3 is here redundant given 
that the simulated gas and galaxies are all at the same epoch 
already. 



5.2.1 Coherent motion of gas and galaxies 

In Fig. [To] (a) we show the Lya-galaxy ^(s) for the Lya in z- 
space and M* ^ 10^ M© galaxies in z-space and real space. 
Fig. [TD] (b) shows the same plot but now for the sample 
of A/t ^ 10* M0 galaxies. Since the high mass sample has 
poorer statistics here we mainly focus on the M, ^ 10* M0 
sample in Fig. [10] (b) but we believe that the same inter- 
pretations apply to both samples. The first thing to note is 
that the anti-correlation is increased as we go from r — 2 to 
2 — 2. This is surprising if we assume that random Gaussian 
motions dominate galaxy peculiar motions. In this model 
we expect more smoothing in 2-space than in real-space and 
that is not observed. 

Fig. [To] (c and d) shows the Galaxy-Lya transmissiv- 
ity correlation function from the Lya in real- and redshift- 
space with (c) M* ^ lO'' M© and (d) M* ^ 10* M© galax- 
ies. Focussing on the low mass sample of M* ^ 10* Mq 
galaxies, the open diamonds show the results from galax- 
ies in redshift-space and the Lya in real space while the 
pink asterisks illustrate the results from galaxies and the 
Lya both in real space. At separations s < 5 h~^ Mpc, we 
see that Galaxy-Lya transmissivity correlation functions in 
z — z (light blue squares) and r — r pink asterisks both lie 
lower than the Lya-galaxy r — z version (open diamonds). 
This behaviour is again different from the expected result 
for random Gaussian motions when the r — z result might 



be expected to lie between r — r and z — z. This and the pre- 
vious anomaly can be explained by coherent motion when 
the galaxies and the gas move together and make the corre- 
lation from the same space (2 — 2 or r — r space) stronger 
than the correlation from r—z space combinations. Although 
the positions of the clouds and the galaxies move, they move 
coherently so that their relative positions in r—r and z — z re- 
main the same, similar to the effect seen in Fig.[T] This would 
leave the cross-correlation function the same when redshift 
or real space was used consistently but undergo smoothing 
when gas was used in z-space and galaxies in redshift-space 
or vice versa. Similar but noisier effects are seen in the sam- 
ple of M, J5 10^ M0 galaxies in the left-hand panels. 

Ignoring the real-z-space combinations, we see that both 
the real-real and z — z results show the same trends. At a 
distance s > 5 /i^^Mpc, the measured (f/fcon) increases 
to reach the mean value. As separations decrease below 
5 h^^Mpc, the transmissivity decreases and hence the Hi 
density increases as we approach the galaxy. 



5.2.2 Dynamical Infall in ^(s) and ^(cr, n) 

As discussed earlier, we note that the mean transmissivities 
for the simulated real and redshift-space Lya forests are 
somewhat different, with the redshift-space forest having a 
mean of T = 0.72 and the real space forest a mean of T = 
0.69. This is surmised to be the result of Lya absorption lines 
overlapping increasingly once redshift-space distortions are 
applied and potentially indicative of infall effects. We shall 
now discuss dynamical infall of the gas further, by looking 
for the effect of gas infall on the cross-correlation function. 

To better visualise any distortions in the cross- 
correlation, we calculate the function ^ = 1 — T/T, using 
a value of T = 0.69 in the real-space case and T = 0.72 
in the redshift-space case. The results for ^(s) are shown in 
panels (e) and (f ) of Fig. \W\ 

We fit the real-space cross-correlation function ^(r) by 
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Figure 10. Lya-galaxy transmissivity, {T(s)) from GIMIC Af* ^ lO'' Mq and Af* ^ 10* Mq galaxy samples. The error bars are 



calculated by using jack-knife method, (a) Af* ^ 10 M© galaxy sample with z— space galaxies 



-I- z- 



-space gas (light blue squares) 



and real-space galaxies -|- z— space gas (purple circles), (b) Same as (a) for Af* ^ 10* Mq galaxies, (c) Lya-galaxy transmissivity for 
Mi, ^ 10® Mq sample with Lya-galaxies in rcal-z (blue diamonds), real-real (pink asterisks) and in z — z (light blue squares). A fit to 
the real-real f(r) with rg = 0.30 h~^ Mpc, 7 = 0.90 is shown (red line). The expected Lyo-galaxy ^{s) in z — z space with I^Lya = 1-3 
>/3gai = 0.35 and {wV^ = -\/95^"-i-T20^ kms~^ is also shown (light blue dashed line). Note that these velocities arc calculated from 
pairs within s < 1 /i~^Mpc. (d) Same as (c) for the A/* ^ 10* Mq galaxy sample. Here the ^(r) fit assumes tq = 0.21 h~^ Mpc, 7 = 0.95 
and the f (s) model assumes /S^^y^ = 1-3 ,/3ga; = 0-53 and {w'^) = \/71^~+T20^ kms~^. (e) ^(s) and ^{s)/^{r) for the high mass galaxy 
sample. The blue line in the lower panel represents the linear theory prediction with fBg^i = 0.35, f^Lya = 1.3 giving i{s)/S,{r) = 1.64. 
(f) Same as (e) for the A/* > 10* Mq galaxy sample with Pg^i = 0.53 , PLya = 1-3 giving C(s)/C(r) = 1.75. 



using C(r-) = (7^) with ro = 0.30 ft"^ Mpc, 7 = 0.90 for 

Mi, > 10^ Mq galaxies and ro = 0.21 h^'^ Mpc, 7 = 0.95 
for M* ^ 10* Mq galaxies. Taking into account a Kaiser 
boost of f(s)/^(r) « 1.5, this ro for the high mass sam- 
ple is reasonably consistent with the sq =0.3 h~^ Mpc fit- 
ted to the GIMIC cr oss-correlation ^(s) in Fig. [10] (see also 
ICrighton et al.ll201ll ). Thus the pink lines in Fig. \W\ (c and 



d) represent T{r) — T — ( — ) ''. T is the mean Lya trans- 
missivity at 2 ~ 3 (A05). 

We then apply t he lin ear theory redshift-space distor- 
tion model of lKaiseri (|l987l ). in the form for cross-correlation 
e(s)/e(r) = (1 + ^(P^ai + I3lv „) + ^Pgai^LyJ discussed by 
iMountrichas fc Shanks! (|2007l ) . As indicated earlier we shall 
assume the value Phya = 1-3 to represen t the 1.0 < 1 3 < 1. 6 
range suggested from the simulations of [McDonald] (|2003h . 
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Th e z-space distortion ana lysis is then implemented follow- 
ing |da^Aiigelaet_alJ (|2008r ). Focussing again on the M* > 
10 M© galaxies sample, we assume the appropriate values 
above of ro and 7, with dynamical infall values Pgai = 0.53 
and PLya = 1-3, predicting ^(s) — 1.75 x £,{r), plus veloc- 
ity dispersion sf{173/V^p'+120^ = ITl kms~\ This pro- 
vides a poor fit even when we reduce the galaxy velocity 
dispersion from 122 kms^^ to 71 kms^^, previously sug- 
gested to be the 'within-clump' value (see Fig. \TU\ (d, f)). 
It is the pLycyi = 1-3 value that is proving too high here, re- 
sulting in this model tending to over-predict the Lya-galaxy 
cross-correlation function. The M* ^ 10* Mq galaxy sam- 
ple rather indicates in (f) that ^(s) « 1.4 x ^(r) implying 
PLya ~ 0.7. Similar results are found for the high mass sam- 
ple of M, ^ lO'"* Mq galaxies (see Fig. \W\ (c,e)) where fits 
in (e) indicate /3lj/ci ~ 0.15. 

The GIMIC Lya-LBG ^{a, tt) results are presented in 
Fig. lllf a.b.c.dV Panels (a) and (b) show the Lya-LBG result 
from M* > 10^ M© galaxies in r — r and z — z respectively 
while panels (c) and (d) present the Lya-LBG result from 
M* ^ 10* Mq galaxies. Panels (e) and (f) give the (w?)^''^ 
versus /3 Ax^ contours for GIMIC Lya-LBG C(<t, tt) with 
Mi, ^ 10® Mq galaxies and panels (g) and (h) the same 
for the Mi, ^ 10* Mq galaxy sample. For the models we 
have assumed the same values of ro, 7 and (in the z — z 
case) Pgai as above. We again note that we see little change 
in the shape of ^(cr, vr) between the r — r and z — z results, 
suggesting that the effect of the peculiar motions, apart from 
coherent bulk motion, is again small. Restricting the fitting 
range to a, tt < 12 h~^Mpc we find a fit to the r — r ^{a, tv) 

giving PLya = 0.02 and {w^) ~ 20 kms^^, with reduced 
X^ = 1.2. Fitting to the z — z correlation function gives 
Phya = 0.15 and (10^) = 40 kms~^, with reduced x^ = 
1.8. Panels (g) and (h) show equivalent contour plots to 
panels (e) and (f) but with M* ^ 10* M© galaxies. The 
best fit for r — r gives PLya = 0.55, {w^) =60 kms~^, 
and reduced x^ ~ 0.6, assuming ro = 0.21 h^^Mpc and 7 
= 0.95. The best fit ior z - z gives PLya = 0.7, (to^)^''^ = 
60 kms~^, and reduced x^ ~ 0.8. 

For the M* ^ 10* M© galaxies, the values that we would 
expect are /3 — (ifz) =0 for the r — r case and PLya = 1.3 

and (lUz)^''^ ^ V1202 + 71^ = 139 kms"^ for the z - z 
case. From the Ax^ contours in Fig. [11] we see that in the 
z — z cases such values are quite clearly rejected by the 
simulated data, the more so when it is realised that the best 
fit parameters for both mass ranges are the same m. z — z 
and in r — r. We conclude again that the main result is 
the lack of difference between the real and redshift space 
fits for both high and low mass galaxies. We conclude that 
the effect of peculiar motions on the simulated z — z 2-D 
correlation function is smaller than expected both for the gas 
infall parameter and for the gas-galaxy velocity dispersion. 
The explanation may be due the gas motion being coherent 
with the galaxies and so (w?) ~ 120kms^^ may need 
reducing to ("w^) ~ 71kms^^ as for the galaxies implying 

(ui^) ~ lOOkms^^ which is less easy to reject at least in 
the case of the low mass galaxies in (h). 



5.2.3 Simulation and observation compared 

We next compare the simulated results for Lya-galaxy ^(s) 
with the Keck-|-VLRS data as shown in Fig. \8\ The Keck 
-I- VLRS result shown is the full sample including A03, 
A05 and VLRS rather than VLRS plus the 6 re-reduced 
Keck QSOs because the larger number of QSOs means that 
the errors are smaller at s < 1 h~^Mpc. The pink dot- 
dashed line shows the GIMIC z — z space result. We have 
assumed the same model as Crighton et al to fit the GIMIC 
resuh ie {r(s)) = T - (s/so)"^ with so = 0.3 /i~^Mpc 
with a small-scale cutoff with T = 0.29 for s < 0.5 h^^ 
Mpc. The blue solid line is the expected Lya-Galaxy ^(s) 
when convolving the empirically determined velocity errors 
of (wl)^'"^ = ^(4202 - 2 X 1002)/2 = 280 kms'^ into this. 
Our bin size is 0.5 /i^^Mpc. The GIMIC model convolved 
with the empirical velocity dispersion lies above the LBG- 
Lya data but only at the \ — la level and we regard it as 
being in excellent agreement with the combined LBG-Lya 
data. 

We finally compare the simulated results for Lya-galaxy 
^((7, tt) now with the 6 Keck QSOs 4- VLRS data as shown 
in Fig. [9l We have seen that the observed best fit pa- 
rameters for the Keck + VLRS data are PLya = 0.02, 
(tu?) — 100 kms^^ and we noted that the measured in- 
fall parameter of PLya = 0.02 was ~ 3(t lower than the 
expected value of PLya = 1.3. The best fit parameters for 
the simulated M* ^ 10* Mq galaxy sample are PLya = 0.70, 



(-.^> 



2\l/2 



60 kms . This simulated PLya value therefore 



lies between the expected value of PLya = 1.3 and the ob- 
served Keck+VLRS value of PLya = 0.02, an intermediate 
result where the Ax^ contours shown in panel b of Fig. [2] 
allow consistency with both these values. We conclude the 
Keck -f VLRS data also prefers low values of PLya like the 
simulation. We also note that neither the minimum value 
of (w^) = 297 kms~^ from LBG velocity error nor the 

(10?) = 320 kms~^ value from including the full simu- 
lated (w^) ~ 120 kms~^ is consistent with the data at 
the few sigma level. We conclude that while the Keck -f VLRS 
PLya and (i«z) estimates are low compared to initial ex- 
pectation from theory, they may still be consistent with the 
similarly low values of these parameters estimated from ^(s) 
and ^((T, tt) in the simulations. 



5.3 Lya auto-correlation 

We measure the Lya auto-correlation along and across QSO 
lines-of-sight. By doing this, we can measure the veloc- 
ity dispersion of Hi gas to see if this has a significant ef- 
fect on the Lya-LBG cross-correlation function. Following 
ICrighton et al.l (|201ll ). for each pixel in a QSO line-of-sight, 
we calculate 



T 

J = i-l 

T 



(7) 



where T and T are the measure and the mean transmissivity. 
We used this to calculate the correlation 



C(Ar) = {5(r)(5(r -f Ar)) 



(8) 
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Figure 11. The top row shows the GIMIC Lya-galaxy 5((t, tt) from M* ^ 10^ M0 (a,b) and M* ^ 10* M0 galaxy samples (c,d) and 
in real-real space (a,c) and in z — z space (b,d). The errors were calculated by using the jack-knife method. The bottom row shows the 
pairwise velocity dispersion ((lo^) )-infall parameter(/3iyQ) Ax^ contours for GIMIC Lya-LBG g((T, tt) fits in the range s < 12 h~^Mpc. 



(e) The best fit for real-real M* > 10^ M0 sample is ^^^^ = 0.02, (w^) 



2\l/2 



20 kms ^, assuming rg = 0.3 h ^Mpc and -y = 0.9. (f) 



The best fit for z - z J\/* > 10^ M@ sample is ^Lya = 0.15, (w^) 



2\l/2 



40 kms~i. (g) The best fit for real-real Af* > 10* Mq sample is 



PLya = 0.55, («)?)i/2 = 60 kms' 
0.70, {w^y^^ = 60 kms-i. 



assuming rg = 0.21 h -"^Mpc and 7 = 0.95. The best fit for z — z M* ^ 10* Mq sample is fiz^yo 



We sum all pixels with the separations Ar, both parallel and 
perpendicular to r. We repeat this process for the QSO data 
listed in Table 1 and for the simulations. 

Fig. [12] (a) shows the auto-correlation of Lya pixels 
along the line-of-sight from the observational data. Keck, 
VLRS and combined samples are presented by pink dia- 
monds, hght blue asterisks and blue circles respectively. Er- 
ror bars were estimated by using jack-knife method. W e first 
compare these to the result from lCrighton et al.ll201ll (trian- 
gles) who measured the auto-correlation using 7 high resolu- 
tion QSOs (resolution FWHM ~ 7 kms~^). The agreement 
appears to be excellent. They all show similar results at small 
scales. We also compared our results with the recent BOSS 
result of lSlosar et al.ll201ll (squares) but it is mostly seen at 
large scales as they are more interested in large-scale corre- 
lations. 

The auto-correlation of simulated Lya pixels are pre- 
sented in Fig. [T2] (b) . The real and redshift-space Lya auto- 
correlations are shown as the asterisks and diamonds, re- 
spectively. The real-space gas simulations are derived from 
the redshift-space simulations by switching off the effect of 
peculiar velocities in the output data. We then compare 
these to the result from Keck-fVLRS. They show a good 
agreement. 

Focussing on the simulation, we see again that the red- 
shift and real-space correlation functions are quite simi- 
lar in amplitude and form. At small scales, convolving the 
real-space correlation function with a Gaussian of width 
120 X \/2 = 170 kms~^ representing the simulation gas 



peculiar velocity (see Fig. ^ is seen to overestimate the 
small-scale turnover in the redshift-space correlation func- 
tion. This is better fitted by a velocity dispersion of width 
30 kms~^ but this assumes I^Lya ~ which is unphysi- 
cal. Now at larger scales we see that the observed ratio of 
^s/£,r ~ 1.5 is smaller than that expected from the linear 
theory in Equation O For the gas the value of the bias is 
b « 0.3 which arises from comparing the mass linear corre- 
lation function to the simulated real-space correlation func- 
tion in Fig. [12] in the range 1 < r < 6 h~^Mpc. This bias 
corresponds to PLya ~ 3.3 which impli es ^{s)/^{r) ~ 5.4. 
But again as noted bv iMcDonaldl (|2003l ). ^Lya has no sim- 
ple relation to density bias as for galaxies. ^Lya has to be 
estimated from simulations and the simulations of McDon- 
ald et al implied a range ^Lya = 1 — 1.6. If we therefore take 
PLya ~ 1.3, then this predicts ^{s)/£,{r) = 2.2 from Equa- 
tion [3] whereas the simulated value in Fig. [12] is more like 
^(s)/^(r) = 1.5 which corresponds to /3 ~ 0.7. The reason 
for this discrepancy is unclear. With ^Lyct ~ 0.7, the best fit 
velocity dispersion is (w^) =70 kms^^. Models where we 

fixed l3Lya = 1.3 and took (iCz) = 170 kms~^ as expected 
from Fig. [5] were also strongly rejected. With PLya = 1.3, 
again a best fit value of (w^) =70 kms~^ was found 
although the model was still rejected in a chi-square test. 

Whatever value of PLya is chosen it appears that 
the fitted value of the velocity dispersion is much lower 
than measured in th e simulation. However, as shown by 
ICrighton et al.l (|201l[ ). the intrinsic width of the Lya lines 
convolved with the instrumental response of the spectro- 
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Figure 12. (a) The auto-corrclation of Lya pixels along the 
line-of-sight. Keck, VLRS and combined samples arc shown by 
pink diamonds, light blue asterisks and blue circles respectively. 
The measurement of Crighton ct al. (2011) is also shown (tri- 
angles). The recent BOSS result of Slosar et al. (2011) is also 
shown (squares), (b)-top panel ; The auto-correlation functions 
of GIMIC Lya pixels at z = 3.06 in the Oct simulation. Real- 
space (pink asterisks) and z — space (blue diamonds) results are 
shown. Errors are calculated via the jack-knife method. A dou- 
ble power-law fit to the real-space {(r) with roi = 0.0049, 71 = 
0-5, ro2 = 0.054, 72 = 1.1 is also shown (pink line). The brown 
dot-dashed line is the expected result for the Lya ^(s) in z-space 
if we convolve in the velocity dispersion {w'^') = 170 kms^^ 
and Pijyci = 1.3 to the z — space distortion model. The blue solid 
line is a predicted ^(s) with (ui?) =70 kms^-'-, /S^^^ = 0.7 
which appears to fit the simulated data. Bottom panel : GIMIC 
^(s)/5(r') with jack-knife error bars. The dashed line corresponds 
to 5(s)/5(r) = 2.2 as predicted from linear theory with Pi^y^ = 
1.3. 



Figure 13. GIMIG Lya g((T, tt) at 2 = 3.06 in real and z-space. (a) 
Lya real-space autocorrelation function, (b) Lya z-space autocor- 
relation function, (c) Pairwise velocity dispersion ((lo^) )-infall 
paramctcr(/3) Ax^ contours for GIMIC, fitting real-space ^(cr, tt) 
in the range s < 8 ft-^Mpc. The best fit is /3 = 0.02, {wl)^^'^ = 
80 kms~^ and reduced x^ = 1.1, assuming roi = 0.0049 /i~^Mpc 
, ''02 = 0.054 fc~^Mpc ,71 = 0.5, 72 = 1.1 and break point is at 
0.40 h-^Mpc. (d) As (c) for z-space. Here the best fit is /3 = 0.02, 
,1/2 



{wl) ' = 80 kms-l and x^ 
power-law fit for 5('r). 



2.4, assuming the same double 



graph can induce artificial autocorrelations at scales < 
0.7 /i~^Mpc, so this eflect may contribute to the poor fit 
of the peculiar velocity z-distortion model on small scales. 

We note that the z-distortion model for the Ly-a auto- 
and cross-correlation assumes spherical symmetry as we 
move from real-space to redshift-space and the Lya auto- 
correlation function involves summing along and across QSO 
lines of sight which may not be exactly spherically symmet- 
ric. However, we shall see that this explanation cannot apply 
to the Lya ^(a, tt) which we calculate next and which gives 
consistent results. 

For each pixel in the Lya line-of-sight, we next calculate 
the Lya ^{a, tt) by using. 



e(a,^) = 



N{a,n) 



(9) 



where < DT{a, n) > is the number of Lya pairs weighted 
by the normalized transmissivity: T^, for each separations. 
N(cr, tt) is the number of Lya pixels that contributed to each 
pair. 

The Lya ^(a, tt) results at z = 3.06 for the Oa simu- 
lation are shown in Fig. 1131 the left hand side showing the 
result in real space and the right hand side showing the re- 
sult in redshift-space. We see similar results for r — r and 
z~z space. Again there is little evidence of smoothing due to 
Gaussian convolution or flattening due to infall in Lya auto- 
correlation. Thus the best fit model for r — r is (wz) = 
80 kms~^, ^Lya = 0.02. For z — z the best fit parameters are 
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the same, Phya = 0.02 and {wlY = 80 kms~^. The ex- 
pected value of jSLya ~ 1.3 is as strongly rejected for z — z as 
for r — r with Xred ~ H)- We again conclude that the effects 
of infall in the gas in the GIMIC simula tion are much less 
than expected from the previous work of iMcDonaldl (|2003l ) 
with an upper limit of PLya = 0.7 from ^(s) and PLya ~ 
from Lya ^(cr, tt). Given /? = 0.7, a gas velocity dispersion 
of (wz) =70 kms~^ fits the simulated Lya ^(s) and this 
is close to the sub-l/i~^Mpc value of the velocity dispersion 
estimated for simulated galaxies, due to correlated motions. 
Good agreement is observed between the GIMIC simulation 
and the Keck+VLRS gas ^(s). We have not calculated the 
Lya 2-D autocorrelation as the QSO sample does not have 
a high enough sky density to properly probe the on-sky pro- 
jected profile. 



6 DISCUSSION 

A feature of this study is that we want to combine the 
power of the VLRS at large spatial scales with the statis- 
tical power of t he K eck sample of A05 at smaller scales. 
ICrighton et al.l (|201ll ) have included the Keck data in the 
LBG-Lya cross-correlation function by simply using an er- 
ror weighted combination of the Keck and VLRS correlation 
functions. Our aim here was to combine the two surveys for 
2-D, ^((T, tt) correlation function analyses at the deeper level 
of the Lya fiuxes and LBG positions. We therefore included 
940 2.67 < 2 < 3.25 LBGs from the Isteidel et af] (|2003l ) 
Keck samples. We also re-reduced 6 high resolution spectra 
of the QSOs in these fields from the ESO or Keck archives. 
With « 3000 galaxies the combined VLRS and Keck surveys 
covering the widest wide range of spatial scales are ideal to 
study the dynamical relationship between galaxies and the 
IGM at z^3. 

The other new feature of this work is our use of the 
GIMIC simulations to help interpret the above correlation 
function results. The simulations are used to create synthetic 
Lya spectra and galaxies. We study both galaxy clustering 
and the relationship between gas and galaxies via the auto- 
and cross-correlation functions in both 1-D and 2-D. 

We first compared the simulated galaxy-galaxy results 
in real- and redshift-space. The LBG-LBG S,{r) and ^(s) in 
real and redshift space appear to be power laws at large 
scales. They show strong clustering at small separations. 
At smaU distances {s < 1 h~^ Mpc), the LBG-LBG C(s) 
in redshift-space tends to have lower clustering than ^(r) 
in real space, while at larger scales the LBG-LBG ^(s) re- 
sults appear to have higher clustering. Qualitatively this is 
as expected from 'finger-of-God' effects at sub 1 h~^Mpc 
scales and dynamical infall at larger scales. Quantitatively, 
the large scale Kaiser boost for the galaxies is consistent 
with expectations. However, at smaller scales the peculiar 
velocity dispersion measured in the simulation overestimates 
the difference between real and z-space correlation functions . 
Similar results may have been found bv lTaruva et al.l (|2010l ) 
who found that at high redshift fitting finger-of-god damping 
terms, as we do here, tended to underestimate the peculiar 
velocity dispersion predicted by linear theory. Certainly, a 
'local' velocity dispersion measured relative to galaxy pairs 
with separations < 1 h~^Mpc produces improved agree- 
ment. For the galaxy ^{a, n) there appears again to be little 



difference between the real-and z-space correlation functions 
but the results are still noisy. Clearly none of the correlation 
function statistics will pick up the most obvious indication of 
peculiar velocity in these simulated data which is the bulk 
motion of « 100 kms^^ across the simulated volume. We 
will see that these unexpectedly small effects of peculiar ve- 
locity are also apparent as we move to study the gaseous 
component. 

We then checked the LBG-Lya transmissivity cross- 
correlation function as calculated directly from spectra and 
LBG positions of the Keck survey against the results of A03 
and found good agreement, although we only had access to 
6 out of their 8 QSO spectra. We checked the results for 
the combined K eck and VLRS s urvey s against the weighted 
combinations of lCrighton et al.l (|201ll ) and again found good 
agreement. Our overall results for the Keck-|-VLRS LBG- 
Lya ^(s) resemble more those of A05 and lRakic et al.l (|2012f ) 
than those of A03 ie a continuous decrease in Lya transmis- 
sivity around the LBG with no evide nce for a spike in trans- 
missivity due to feedback being seen. ICrighton et al.l (|201ll ) 
noted that such a spike could still be present but smoothed 
away by the errors in the LBG velocities. 

The inclusion of the gas component along with the 
galaxies makes possible the search for gaseous infall. As- 
suming the appropriate dynamical infall parameter for the 
simulated M* ^ 10* M© galaxies, the gas- galaxy cross- 
correlation model gave PLya ~ 0.7, lower than the expected 
range of 1 < ^Ly^ < 1.6 from iMcDonaldl (120031 ). This small 
effect of the peculiar velocities was confirmed in the analy- 
sis of galaxy-Lya ^{a, n) where again the circular contours 
of the real-real ^{a, n) was almost matched by those from 
the z-space version. Again there seems consistency with our 
view of a broadly coherent motion between gas and galaxies. 

From the Lya auto-correlations ^(s) and ^(a, tt), we 
see similar results with again small differences between 
real- and redshift-space. At small scales, the velocity dis- 
persion needed to fit the simulated ^(s) is less than mea- 
sured directly in the simulation, although this may be 
partly explained by the intrinsic width of the Lya lines 
contributing artificial autocorrelation below separations of 
< 0.7 /i~^Mpc. At larger scales, the value of C(s)/^(r) is 
more co nsistent with PLya ~ 0.7±0.3 rather than the range 
given bv lMcDonaidI (|2003l ). PLyc ~1- 1.6. 

At larger scales, one possibility to explain the low gas 
infall rate may be due to the presence of feedback in the 
GIMIC simulations. Galaxy- wide winds of 600 kms~^ are 
invoked in the GIMIC simulations and this is a significant 
amount since this corresponds to 6 /i~^cMpc. These winds 
are modelled by each star particle that forms, imparting a 
randomly directed 600 kms^^ kick to 4 of its gas particle 
neighbours. It is possible that this outfiow of the gas could 
cancel out some of the expected gravitational infall particu- 
larly in the neighbourhood of a galaxy. However, it remains 
to be seen whether enough gas particles are outfiowing to 
explain the lack of infall in the gas cross- or auto-correlation 
functions. If the effects of gas outflow were detectable in the 
gas dynamics this could be a powerful probe, since there is 
not much evidence of feedback from any spike in transmis- 
sion due to lower neutral gas density close to the galaxy. 
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7 CONCLUSIONS 

We have analysed the interaction between IGM and galaxies 
at 2 ~ 3 using VLRS data and KECK LBGs, also the spec- 
troscopy for 17 QSOs in the redshift range 2.5 < z < 4.0 
obtained from the publicly available QSO spectra and a 
spectrum from our own observations using the X-Shooter 
instrument on the VLT (Q2359-I-0653). Apart from the ob- 
servational data, we employ the SPH GIMIC simulations to 
analyse the clustering of gas and galaxies. 

1. The observed galaxy real space autocorrelation func- 
tion is more consistent with that measured for simulated 
Mi, > 10^ Mq galaxies than low mass galaxies. When an 
empirically determined velocity dispersion dominated by ve- 
locity errors is convolved with the simulated real space cor- 
relation functions, a similar preference is found for the lower 
amplitude clustering of the M* ^ 10® M© galaxies. In the 
simulated data the difference between the real and redshift- 
space correlation functions is too small to be self-consistently 
explained by the measured peculiar velocity distribution. We 
suggest that this is the consequence of a scale dependence 
in the measurement of the peculiar motions and that the 
peculiar motions taken within < 1 h~^Mpc give a more con- 
sistent result. 

2. We have checked for the existence of the transmis- 
sion spike near star-forming galaxies in the data and GIMIC 
simulations which would be indicative of the effects of star- 
formation feedback on the IGM. For the data, we combined 
the full VLRS and Keck LBG-Lya datasets to study both 
^(s) and ^((T, tt) and the LBG-Lya correlation functions. 
No detection of feedback is seen in the galaxy-gas corre- 
lation function in the VLRS-|-Keck data. We find that the 
gcis transmissivity monotonically drops towards the galaxy, 
consistent with the density of neutral gas rising towards the 
galaxy position. Although the simulation transmission rises 
when LBG velocity errors are taken into account, the sim- 
ulated and observational results remain in good statistical 
agreement. 

3. The redshift-space gas-galaxy cross-correlation func- 
tion in the simulation is close to the real-space correlation 
function and to some extent this is expected from linear 
theory applied to the Lya forest flux which has a non-linear 
relation with optical depth and thus implies lower rates of 
dynamical infall of gas into galaxies than would otherwise 
apply. We have also considered whether galaxy-wide out- 
flows may be cancelling out the infall effect. 

4. The observed Lya autocorrelation function is also 
consistent with the simulation. At small scales the difference 
between real and redshift-space correlation functions in the 
simulation is again less than expected given the peculiar 
velocity distribution but may be consistent with previous 
simulations of the Lya forest. At larger scales the effects of 
dynamical infall are in line with linear theory, if the non- 
linear relation between flux and optical depth is taken into 
account. There may also be some residual effect from gas 
outflows cancelling out the effects of dynamical infall. 

5. In the simulations, both gas and galaxies show evi- 
dence of a strong bulk motion. This bulk motion is unde- 
tectable by observable correlation functions but may have a 
connection with the local coherence needed to explain why 
distribution of peculiar velocities overestimates the flnger- 
of-God effect. 
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